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INTRODUCTION 

Since the landmark paper by Bellman and Zadeh in 1970 [1], fuzzy sets have been 
used to solve a variety of optimization problems. In mechanical and structural design, 
however, relatively few applications of fuzzy optimization have appeared in the literature, 
even though the formulations and algorithms developed in management science are readily 
applicable in engineering fields. In engineering design fuzzy optimization has most often 
been suggested as a mechanism to represent 'soft' constraints (e.g. [2]) and as an 
alternative formulation to solve multiobjective optimization problems (e.g. [3,4].) A 
formulation for multilevel optimization with fuzzy objective functions is presented in this 
paper. 

With few exceptions, formulations for fuzzy optimization have dealt with a one- level 
problem in which the objective is the membership function of a fuzzy set formed by the fuzzy 
intersection of other sets. This concept dates back to the original work by Bellman and 
Zadeh in which objective and constraints are handled in identical fashion: if f(x) is the 

function to be minimized under constraints g; (x)<0, i=1 ,2 m, the fuzzy optimization 

problem maximizes the membership function po of 

G=GonGin,..,nGm with pG=min(po,pi Pm) 

The set Go has membership function po such that 

po(f)-» 1 as f decreases 
while, for other Gj's, 

Pj (gj ) 1 as gj decreases from 0, and 

Pj (gj ) — > 0 as gj increases from 0. 

This model has been used extensively to relax ('fuzzify') the constraint set and to deal with 
multiple objectives. 

A somewhat different problem is discussed here. First, the goal set G is defined in a 
more general way, using an aggregation operator H that allows arbitrary combinations of 
set operations (union, intersection, addition) on the individual sets Gj. This is a straight- 
forward extension of the standard form, but one that makes possible the modeling of 
interesting evaluation strategies. This feature has been discussed in detail elsewhere [3,4] 
but, for completeness, it will be briefly outlined in the next section. 

A second, more important departure from the standard form will be the construction of 
a multilevel problem analogous to the design decomposition problem in optimization [5-8]. 
This arrangement facilitates the simulation of a system design process in which 

- Different components of the system are designed by different teams. 

- Different levels of design detail become relevant at different time stages in the 
process: global design features early, local features later in the process. 
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SINGLE -LEVEL PROBLEM 

The optimization problem can be solved in a single level when a design alternative 
can be fully described and evaluated using local design variables and functions. Let xe R n 
represent the vector of design variables. The goal of the problem is to find the design x 
that 


maximizes PG=h(pi,P2--.Pp), 
the membership in the design goal 

G=H(Gi,G 2 ,..,G p ), 
subject to crisp constraints 

xe X= {xe R n : Q\ (x)<0, i=1 ,2 m}. 

Each individual goal Gj has membership \i\ and represents a desirable design quality such 
as strength, low cost, or reliability. The mapping from x to each of the pfs depends on 
specific details of the design process and designer's preferences and it is assumed to be 

known. The function h:[0,1]P-»[0,l] is an acceptable connective associated with H, as 
defined in [3,9]. For instance, in a standard 3- objective fuzzy optimization problem 

H(G 1 ,G 2> G 3 )=GinG 2 nG 3 and h(pi,p 2 ,p 3 )=min{pi,p 2 ,p 3 }. 

This is by no means the only possible or meaningful choice for H and h. For instance, the 
choice 


H(Gi,G 2 ,G 3 )= (GinG 2 )u(GinG 3 )u(G 3 nG 2 ) 
indicates that the optimum solution needs to satisfy only 2 out of 3 objectives. The 
connective h can be built using different intersection and union operators. An extensive 
review of these connectors can be found in [10]. 





0=6,063 n G 3 G^G^GjIulGnGgtufGnGj) 
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MULTILEVEL PROBLEM 

In a multilevel approach the system being designed is decomposed into several sub- 
systems. This decomposition is often made along boundaries determined by the 
organization and expertise of design teams involved in the design task. Each team becomes 
responsible for the design of its own sub- system and is an expert only in this limited area. 
Sobieszczanski-Sobieski, J. [5,7,8], Haftka [6], and Parkinson et. al. [1 1], among others, 
have suggested different approaches to the optimization problem in this setting. 

Two important questions arise in the multilevel optimization problem: what is a 
meaningful objective function to guide the optimization and how to maintain feasibility when 
only local constraints can be evaluated exactly. The fuzzy set approach offers an 
operationally useful answer to the first question. The second question will be addressed 
using convex approximations to non- local functions. 

In the fuzzy multilevel problem each sub- problem (j) has its own variables x(D and 
goal GO) . The sub- problems contribute their goals to the overall design objective G, formed 
in a hierarchical fashion using set operations on the sub-problem goals. This arrangement: 

- Introduces the flexibility of the fuzzy formulation at the sub- problem level. All of the 
tools from fuzzy optimization are available there, including the traditional 'crisp' optimization 
problem. 

- Facilitates the comparison of dissimilar objectives arising in different environments: 
membership functions are dimensionless. 

- Makes possible the construction of a design goal that varies as the design process 
unfolds. The vagueness that characterizes the early stages of the design can be easily 
modeled at the local level in the definition of the local goal. 








GLOBAL PROBLEM 

Let xe R n be the vector of all design variables. The solution to the global optimization 
problem 

maximizes piG(x)=hG(|x( 1 ) ,|o.( 2 > 

under the constraints of the local problems defined below. The functions p(i) are 
membership functions of the sets G(i) corresponding to acceptable solutions to thej-th 
sub- problem and the function p G corresponds to the fuzzy set G that describes the 
optimum global solution, 

G=H g (G(1),G( 2 ),...,G(p)) 

The operators h G and H G describe the way different sub- problems interact and, in general, 
change to reflect different stages in the design problem. 

LOCAL PROBLEM PG) 

Let x(i)e R n (i) represent the vector of local design variables and let y(i)e R N n (i) 
represent design variables outside the j-th problem, i.e., yG)=x- x(i). Problems are assumed 
to be coupled: objectives and constraints in the j-th problem are functions of y(i). However, 
while x(i) varies in the j-th problem, y(i) remains fixed. The goal of the local problem is to 
find the design x(i) that 

maximizes p G (x(i);y(i),f(j)(x(i),y(i))) 
subject to crisp constraints 

xG)e X(i)= {x(i)e R n (i) : g(i)j (xG)yG),fG)(xG),yG)))^0, i=1,2 m(j)} 

Local goals are coupled with other sub-problems via y(J’) and f(i). The fG'Vs are functions 
evaluated outside problem PG). Convex approximations are used to evaluate these 
functions. 
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EVALUATING NON-LOCAL FUNCTIONS: Maintaining feasibility. 

A challenging difficulty in the multilevel problem involves maintaining global 
feasibility while solving the sub-problems. The question is, how to insure that aH variables, 
not just the local variables, remain feasible without actually evaluating non-local 
constraints. An answer to this question will be attempted using convex approximations of 
non-local functions. 

Convex approximations were introduced by Starnes and Haftka [12] and used by 
Fleury and Braibant [1 3] to solve a range of structural optimization problems, lff(x) is a 
differentiable function in R N , the convex approximation of f at x° is the function 



where the (+) sum is over positive derivatives while the (-) sum is over the negative 
derivatives of f. It is easy to show that for positive variables the convex approximation is the 
more 'conservative' approximation out of the linear, reciprocal, or concave approximations of 

f [12] . Indeed, the set f(x)<0 is often contained in the set * ( x ) - ® . This motivates the 
following variation of the problem P(i): 

Find x(i) eR n (i) that maximizes fir; subject to the local constraints 


/ j\ (i) 

c (J) e X 


-{ 


(j) 

x e R :g 


:g(, J) (x (l, ;y (i, y j, (x (i> l y (I) )) <0, i = 1,..,m(j)} 


and the convex approximations of the non-local constraints at the starting point, 


(n ~ (k) r 
x e X = { 


x (J) e R (n<J)> :g' k> (x (J) ) <0, i= 1,..„m(k)},k = 1 p;k* j 


Along with a step-size restriction, the addition of these constraints is often enough to insure 
that x(J) remains within all feasible sets, local and non-local, without need to evaluate the 
non-local constraints g( k ) and functions f (i) exactly. The convex approximations are easy to 
compute and only require the kind of sensitivity information that is often available in many 
nonlinear programming algorithms. Efficient methods to compute these derivatives have 
been reported in the literature [14]. 



TRUE CONSTRAINT 
SET 
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SEQUENCE OF SOLUTION 


In the solution strategy of the multilevel problem one local problem, the active problem 
P(i), is solved at a time by searching only in the space of local variables xO). The choice of 
operator he determines the relative importance of each sub-problem and hence the order in 
which problems will be solved. Consider for example a non-associative connective such as 

h G (p(i),p( 2 ))=min{p(i),p(2)} ( G=G( 1 )nG( 2 )) 

If pC >)>p( 2 ) at a given point x, the global problem is dominated by the sub-problem goal 
G( 2 ) at that point and efforts should be directed to improve p( 2 ). If, on the other hand, 
hG (|x( 1 ),|i( 2 ))=p.( 1 ) x p( 2 ) 

some gain is possible even if only pO) is improved. Indeed, sub-problem goal GO) 
dominates with possibility p( 2 ). In general, goal GO) will dominate with possibility 


(+)=right, (-)=left derivatives 
To simplify the solution strategy only the 'min' operator will be used to connect sub- 
goals (any operator can be used within each problem). With this simplification mj is always 
1 for the dominant goals (0 otherwise) and the procedure is simplified. 

In some steps, p(J) can be improved within XO). It may happen, however, that an 
increase in a dominant pO) is possible only after relaxing some constraint that depends on a 
non-local variable y(i). When this happens it may be necessary to make another problem, 
say P(k), k*j, active and seek to relax the troublesome constraint in X( k ) instead of improving 
its natural local objective. A robust strategy to select the next active problem and its objective 
function is essential to the success of the approach. Efforts are being directed toward this 
goal but more research is still required. An example that applies to the 'min' operator is 
outlined below. 


m . = max 

J 


fdh Q \ 

7 $) 


fd h, 

l 


<j) 


At x°, after solving P', M={j: pj=1} S)g= feasible descent direction for pg in Xj 
SJ = feasible descent direction for pi in Xj 

j MAX p* 

IF S G *0, je M, THEN x j 


j MAX p j , 

IF S Q * 0, any j, THEN J 


IF S *0,j e Ml j *i THEN Set relaxation: r k>0» keM 

r k=0 > otherwise. 

MAX MIN{n ,k) + r k } 

X J k 

ELSE for j e M, k‘ = Tk s M: g k (x) = 0, (vV, V k g k ) -> mini 

L kg M J 

MIN g k r , n G (x) >n^x°) 

k 
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EXAMPLE 

The multilevel problem with fuzzy objectives can be illustrated by the following simple 
model of support system for a heavy piece of machinery. The problem is modeled by a 
flexible beam supported by three spring- dashpot components. To illustrate the approach the 
device will be decomposed in two subsystems: the suspension (springs) and the structure 
(beam). 

1 . Suspension: Spring stiffnesses are to be selected to transmit a periodic load to the 
foundation. The amount of force transmitted and the maximum deflection are relevant 
performance measures. In this problem the flexibility of the structure is ignored. 

1.1 At a detailed level , the springs themselves are to be designed. Stiffness, stress, 
clearance, and natural frequency are relevant measures. Coil and wire diameters are 
design variables. 


2. Structure: Beam cross-sectional properties are to be selected to support a load. 
Amount of material, stresses and deflections are relevant performance measures. In this 
problem dynamic effects are ignored. 







The 'min' intersection operator will be used to connect the sub-goals, i.e, the overall 
goal is 

G=Gsusp^Gstmct<^Gspring. PG=min{(i. susp » Instruct. Pspring} 

The objective in the suspension sub-problem (1) is the maximum force transmitted to 
the foundation by the three springs. The weight of the beam is the objective in the structure 
sub-problem. Although these are rather simple objectives for a fuzzy optimization problem, 
no significant loss of generality is introduced. Membership functions for these goals are 
shown below. 

Problem 1.1 is special: it corresponds to the detailed design in the suspension design 
problem.The precise description of the spring (problem 1.1) becomes relevant only after 
some knowledge of the required stiffnesses is available (problem 1). This organization is 
present in typical design problems in which detailed design decisions are taken later in the 
design process. In this problem diameters D and d must be selected so that the spring 
stiffness matches the stiffness k prescribed by the problem (1) without exceeding limits on 
stress. The goal is set as 

Gspring*GstiffnGstres3, Pspring=min{|x stiff . Pstress} 

where 


G s tiff={(D,d): k(D,d)*k} 

Gstress={(D»d): o(D,d)<=OY} 

Membership functions for these sets are shown below. 


k m k a * M'susp 

1* 2 o 

8 n (h, k) < 8max 

Knin- k - ^max 


hT?^h ^struct 

h ,' h 2' h 5 

<*h,k) <<*, 

8(h,k)£8 max 
h_, <h £ h 


maxu . 

D.d r ^ s P nn 9 

freq(D,d) >f min 

^min“ ^ “ ^max 

cl . <d<d, 


min 


max 


min 


'max 



|k(U,d)-k|/k 
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Solution history for the problem is shown in the figure below. 

(0) p susp dominates. Problem 1 is solved with objective p q. 

(1) p susp = Pstruct- Problem 2 is solved. Objective is p s truct and p SU sp is not allowed to 
decrease (G S usp to worsen) more than 20%. 

(2) p susp< Pstruct- Problem 1 i solved with objective p g- 

(3) p spring= Psusp- Problem 1 .1 is solved with objective p spring- At the end of this step 
tighter restrictions on a are imposed ( (*) in membership function above). 

(4) p susp< Pstruct and a non-local constraint (maximum beam deflection) is active 
(figure). Problem 2 is solved to reduce this constraint. Notice that the convex approximation 
is effective to determine constraint activity. 

(5) p susp= Pstruct- Problem 1 is solved with objective p g- 



K 1 [-*-500 #/in] 
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CONCLUSIONS 

The introduction of fuzzy objective functions into the multilevel optimization problem 
adds a new dimension to the problem. Interactions among sub-problems can be 
represented as set operations on fuzzy sets, a feature that introduces flexibility and insight 
into the problem. 

A more robust strategy to select the next active problem and objective needs to be 
developed. Heuristic approaches to determine constraint activity may prove useful here. 

Some of the problems associated with keeping solutions feasible can be handled 
using convex approximations of non-local constraints. This is a promising approach, but one 
that needs further tests. More research is needed on issues such as how often sensitivity 
analyses must be performed to keep approximations accurate and how to determine 
appropriate step size limits. 
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